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ABSTRACT 

Context. It is commonly believed that millisecond radio pulsars have been spun up by transfer of matter and angular momentum from 
a low-mass companion during an X-ray active mass transfer phase. A subclass of low-mass X-ray binaries is that of the accreting 
millisecond X-ray pulsars, transient systems that show periods of X-ray quiescence during which radio emission could switch on. 
Aims. The aim of this work is to search for millisecond pulsations from three accreting millisecond X-ray pulsars, XTE J1751-305, 
XTE J1814-338, and SAX J1808.4-3658, observed during their quiescent X-ray phases at high radio frequencies (5 -H 8 GHz) in 
order to overcome the problem of the free-free absorption due to the matter engulfing the system. A positive result would provide 
definite proof of the recycling model, providing the direct link between the progenitors and their evolutionary products. 
Methods. The data analysis methodology has been chosen on the basis of the precise knowledge of orbital and spin parameters from 
X-ray observations. It is subdivided in three steps: we corrected the time series for the effects of (I) the dispersion due to interstellar 
medium and (II) of the orbital motions, and finally (III) folded modulo the spin period to increase the signal-to-noise ratio. 
Results. No radio signal with spin and orbital characteristics matching those of the X-ray sources has been found in our search, down 
to very low flux density upper limits. 

Conclusions. We analysed several mechanisms that could have prevented the detection of the signal, concluding that the low lumi- 
nosity of the sources and the geometric factor are the most likely reasons for this negative result. 

Key words. Pulsar: general - Pulsar: individual (XTE J1751-305, XTE J1814-338, SAX J1808.4-3658), neutron star, X-ray binary 



1. Introduction 

It is commonly believed that accreting millisecond X-ray pul- 
sars (AMXPs), transient binary systems hosting a fast spinning 
(^spin ~ 1 rns) and weakly magnetised (B ~ 10 8 " 9 Gauss) neu- 
tron star (NS), are the progenitors of the radio millisecond pul- 
sars (MSPs), as argued by the recycling model (Alpar et al. 1982; 
Bhattacharya & van den Heuvel 1991). This model asserts that 
the NS in these systems is spun up by the transfer of matter and 
angular momentum from its low-mass (M < 1 M ) companion 
star via the formation of an accretion disk. When this process 
ends, the NS switches on as a radio MSP. 

A basic requirement for the switching on of the radio emis- 
sion is that the space surrounding the NS has to be free of matter, 
a condition that can be fulfilled during the quiescence phase of 
AMXPs. For this reason this phase constitutes the most promis- 
ing one to investigate for confirming of the recycling model. The 
aim of this work is, in fact, to search for radio millisecond pul- 
sations from a sample of AMXPs in their quiescence phase: a 
positive result would unambiguously establish that AMXPs are 
the progenitors of at least some of the radio MSPs. 

In the past decade, from 1998 April, the date of the discov- 
ery of the first AMXP, SAX J1808.4-3658 (Wijnands & van der 
Klis 1998), several attempts have been made to obtain this con- 
firmation (e.g. Burgay et al. 2003), but, despite the eleven ad- 
ditional systems discovered since (e.g. Galloway et al. 2002; 
Casella et al. 2008), we have not obtained a positive result yet. 



A possible explanation of some of these failures has been 
given by Burderi et al. (2001), asserting that detection of radio 
pulsations from AMXPs could be hampered by matter surround- 
ing the system. In fact, during the so-called radio-ejection phase, 
the pressure of the rotating magneto-dipole could prevent the in- 
falling matter from the companion of the NS in the binary system 
from reaching the NS Roche lobe, forcing it to leave the system 
from the Lagrangian point L\ . Once this happens, such matter, 
which is still carrying the angular momentum, will rotate around 
the two stars embedding the system. Even if the pulsar radio 
emission was switched on and the system in X-ray quiescence, 
it could be absorbed by a free-free mechanism. Since the optical 
depth for the free-free absorption, , depends on the square in- 
verse of the frequency, observations at higher frequencies could 
encompass this problem. 

With this scenario in mind we have undertaken a campaign 
of observation at high radio frequencies (see Table 1) for four 
AMXPs, XTE J0929-314, XTE J1751-305, XTE J1814-338, 
and SAX J1808.4-3658, in their quiescence phase. Results for 
XTE J0929-314 were presented in Iacolina et al. (2009). Here 
we used the same equation to estimate Tff for the other three 
sources, and we obtained the values listed in Table 1. The pa- 
rameters are: for the four sources, m\ ~ 1.4, X = 0.7, Y = 
0.3, y~ 1,T 4 ~ l;forXTEJ1814-338,m_io ~ 1.2, m 2 -0.17, 
P h = 4.27 (Markwardt & Swank 2003); for XTE J175 1-305, 
m_io ~ 8.5, m 2 ~ 0.015, P h = 0.71 (Markwardt et al. 2002); for 
XTE J0929-314, m_ 10 ~ 2.9, m 2 ~ 0.02, P h = 0.73 (Galloway 
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et al. 2002); and for SAX J1808.4-3658, m_ 10 ~ 3.7 (which is 
the average mass transfer rate for the 1998, 2000, and 2002 out- 
bursts, see Papitto et al. 2005), m 2 ~ 0.05, P h = 2.01. In this ta- 
ble, the values of Tg obtained at 1 .4 GHz (the typical frequency 
used to observe pulsar) are much higher than unity, so the ra- 
diation would be totally absorbed, while for frequencies higher 
than 5 GHz we obtained values smaller or close to unity for all 
the sources and, reasonably assuming that the matter enclosing 
the system is clumpy (i.e. there are favourable directions where 
the optical depth is lower than the average values of the Table), 
we have a higher probability of detecting the radio signal. 



Table 1: Optical depth at various radio frequencies, for the four sources. 



Source 




Tg 




1.4 GHz 


6.5 GHz 


8.5 GHz 


XTE J1814-338 


2.2 


1.2 


0.05 


XTE J1751-305 


27.8 


1.1 


0.7 


XTE J0929-314 


5 


0.2 


0.1 




1.4 GHz 


5 GHz 


6 GHz 


SAX J 1808. 4-365 8 


9.3 


0.7 


0.5 



The unknown inclination of the systems negligibly affects 
the estimate of the optical depth, considering that, while the es- 
timated companion mass only slightly increases for a decrease 
in the inclination, the amount of matter along the line of sight 
significantly decreases, because the disk is in the orbital plane. 

2. The sources, the observations, and the data 
analysis method 

Two series of radio observations were taken on 2003 December 
20-23 for XTE J175 1-305 and XTE Jl 814-338, and on 2002 
August 5-7 for SAX J1808.4-3658, using the Parkes 64 m ra- 
dio telescope. Observation parameters are listed in Table 3. The 
data analysis methodology was chosen on the basis of the precise 
knowledge of the spin and orbital parameters from X-ray obser- 
vation. It is the same as was adopted for the source XTE J0929- 
314 presented by Iacolina et al. (2009), where it is described in 
detail. 

The original ephemerides were published by Markwardt 
et al. (2002) for XTE J175 1-305, by Markwardt & Swank 
(2003) for XTE J1814-338 and by Chakrabarty & Morgan 
(1998) for SAX J1808.4-3658, and subsequently refined by 
Papitto et al. (2007, 2008) for XTE J1751-305 and XTE J1814- 
338, and by Hartman et al. (2009) (but also by Burderi 
et al. 2009, Di Salvo et al. 2008 and Hartman et al. 2008) 
for SAX J1808.4-3658. Table 3 reports the most updated 
ephemerides used in this work. 

The first part of the data analysis was to correct the time se- 
ries for the dispersion effects of the ISM; the steps, the ranges of 
DMs used, and the values of the local DMs for the three sources 
are indicated in Table 3 (for SAX J1808.4-3658, we considered 
the highest mass transfer rate assumed in quiescence, m ~ 10~ 9 
M /y, proposed by Di Salvo et al. 2008). We then barycentred 
the data series to correct for the orbital effects considering the 
propagation of the uncertainties in the ephemerides derived from 
X-ray observations over the time range between X-ray and ra- 
dio observations. This time range corresponds to -20000 orbits 
for XTE J1751-305 and -1000 orbits for XTE J1814-338. For 
SAX J1808.4-3658, the adopted X-ray ephemerides are derived 



from the analysis of the secular evolution reported by Hartman 
et al. (2009) (see the Table 3), which refers to the time of -100 
orbits after the radio observations. 

For XTE J1751-305, only the propagation of the orbital pe- 
riod error (90% confidence level) affected the time series losing 
the possibility of detecting the signal: i.e. producing a broaden- 
ing of the pulse of 0.9 in pulsar phase. To reduce this broadening 
to at most 0. 1 in pulsar phase, one has to search for the signal 
at 18 trial values of the orbital period, 9 above and 9 below the 
nominal value, covering the uncertainty range. 

For XTE J1814-338 and SAX J1808.4-3658, the uncer- 
tainty in all the parameters within the 90% confidence level for 
XTE J1814-338 and lcr level for SAX J1808.4-3658 does not 
affect the detectability of the pulsation. We then corrected the 
data series by only using the nominal values of the parameters. 

The last step in this search is to fold the time series using the 
spin parameters reported in Table 3. The spin period range ex- 
plored has to be chosen by considering the nominal value of the 
spin period at the epoch of the radio observations, as explained 
by Iacolina et al. (2009). 

For SAX J1808.4-3658, which is the only AMXP that 
showed more than one outburst, we considered the value of the 
spin period at the epoch of the X-ray observations, Px, derived 
by Hartman et al. (2009) and the value of its derivative, Px, mea- 
sured during the outburst closest to the radio observations time, 
reported by Burderi et al. (2006) and listed in Table 3, which 
turns out to be higher than what is derived by the analysis of the 
secular evolution obtained by Hartman et al. (2009) and consti- 
tutes an upper limit in the search (see below). The folding trial 
values are indicated in Table 3. 

To check the plausibility of the adopted spin period interval, 
for XTE J18 14-338, we derived the spin period derivative, P dip , 
through an estimate of the surface magnetic field, considering 
the magnetic torque acting on the neutron star, as discussed by 
Iacolina et al. (2009). The value obtained is P dip = 6 x 10~ 20 , 
which is lower than Px and then safely contained in our interval 
of spin period trial values. For SAX J1808. 4-3658, Burderi et al. 
(2006) calculated a value of the magnetic field equal to Z?s ~ 
(3.5 ± 0.5) x 10 8 Gauss. The spin period derivative, therefore, 
results Pdip = 5.5 x 10~ 20 , which is once again < P x . 

For XTE J175 1-305 Wijnands et al. (2005) derived an es- 
timate of the magnetic field during its quiescence phase con- 
strained to be < (2.5 - 6) x 10 8 G, using the value reported in 
Table 3 for the distance. The spin period derivative resulted to be 
< (2.7 - 15) x 10~ 20 which is lower than P x , consequently, the 
adopted interval of spin periods is again safely large. 



3. Results 

At the end of the three steps of analysis we produced -67000 
folded profiles for XTE J1751-305, -330 for XTE J1814-338, 
and -200 for SAX J1808.4-3658, reporting the results from the 
folding of the dedispersed, deorbited, and barycentred time se- 
ries. The best ones were displayed for visual inspection to search 
for possible pulsar suspects. 

For XTE J1751-305, the highest S/N obtained was 6.78 cor- 
responding to a peak at Act significance, the corresponding plot 
is shown in Fig. 1 (upper panel), where the grayscale on the left 
shows the signal in the 255 subintegrations in which the whole 
observation was split vs the spin phase, while the parameters 
used for the folding are indicated on the right. The diagram at 
the bottom displays about 4 phases of the integrated pulse pro- 
file. The parameters for the pulsar suspect are DM = 735.41 pc 
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cirT 3 , P olb = 2545.34361 s, v s = 435.317974 Hz. This peak has a 
40% probability of not being randomly generated over the 40755 
trial foldings of the dataset corresponding to one of the two ob- 
servations at 8.5 GHz. Analysing the behaviour of the S/N as a 
function of the spin frequency and DM trial values, a decreasing 
trend from a peak at S/N = 6.78 is found. Since this peak was at 
the limit of our search interval in spin frequency, we considered 
it appropriate to investigate for 11 more steps in the spin fre- 
quency, discovering, in this way, the whole trend: the maximum, 
which corresponded to the same S /N, is defined in all the di- 
rections. The result is displayed in Fig. 1 (lower panel), where in 
the y-axis we have the spin frequency (10 steps above and 1 1 be- 
low the nominal value, corresponding to 0) and in the x-axis, the 
DM (60 steps corresponding to the interval between 485.22 and 
932.54 pc cm" 3 ) plotted for an orbital period P orh = 2545.34361 
s. A clear maximum is well defined, hence supporting the plau- 
sibility of the suspect. 

This result is not confirmed in the other observation (at the 
same frequency) elaborated with the same parameters, but this 
could be due to the dumpiness of the matter around the system. 
This suspect may thus deserve additional investigation in the fu- 
ture. 

The highest S /Ns reached for XTE J1814-338 and for 
SAX J1808.4-3658 were 4.67 and 4.41, corresponding to a peak 
at 2.6<t and 2.2cr significance, respectively, with a probability of 
being real radio signals and not produced by noise of 18% and 
20% over the 186 and 57 trial foldings of their single dataset, 
respectively. Unfortunately, visual inspection of these two re- 
sults and others at lower S /Ns did not provide any evidence of 
the pulsed signal, and no other investigations have shown the 
positive signs in the trend of the S/N obtained for XTE J1751- 
305. Finally, the observations at the other frequency, folded with 
the same parameters, did not displayed any suitable shape for 
the signal. For these two sources, we can thus conclude that no 
pulsed radio signal has been detected at their spin period in our 
observations. 



3. 1 . Upper limit on the flux density 

Considering the parameters indicated in Table 3 and a duty cycle, 
W/Ps, of 15%, we estimated the upper limit on the flux density 
for the three observed sources at the nominal DM with the Eq. 
9 in Iacolina et al. (2009). Including results for XTE J0929-314 
(Iacolina et al. 2009), we obtained the values indicated in Table 
2. 



Table 2: Flux density upper limits for the four sources at the analysed 
frequencies. 



Filename: XRT003 0121 S07 N28%DMR0084-0084.DEO 

Folding Parameters 



Source 


S max 


0"Jy) 


6.5 GHz 


8.5 GHz 


XTE J1814— 338 


52 


25 


XTE J1751-305 




31-30 


XTE J0929-314 


68 


26 




5 GHz 


6 GHz 


SAX J 1808. 4-365 8 


59 


59 



Notes. For XTE J175 1-305 the two observations are both at 8.5 GHz. 
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Epoch of pars (MJD): 
Freq (Hz): 
Freq Der (Hz/s): 
Freq 2ndDer (Hz/s*s): 
Period (ms): 
Period (s/s): 
Period 2ndDer (s/s*s): 



52993.06600000 

435.31797400 

0.00000E+00 

0.000E+00 

2.297171 

0.00000E+00 

0.000E+00 



Observation Parameters 



Epoch at start ofobs(MJD): 
Freq at start ofobs(Hz): 
Period at start ofobs(ms): 

Number of subintegrations: 
Length of subintegration (s) 
Length of data (s): 
Number of bins: 
Length of one bin (ms): 
S/N: 



52993.0659143518 
435.317974000000 
2.297171 
0.0600 
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Fig. 1: Upper panel: Plot with the highest S/N for XTE J1751-305. 
Lower panel: S/N in function of spin frequency and DM. See the text 
for further explanations. 



4. Discussion 

In this section we discuss the results of our findings, including 
results for XTE J0929-314 presented in Iacolina et al. (2009). 
A part for the result obtained for the 8.5 GHz observation of 
XTE J175 1-305, deserving additional investigation, no radio 
pulsed emission has been found above the reported upper lim- 
its in the data analysed. Assuming that the radio emission was 
switched on during the observations of the three sources, we in- 
vestigated the possible reasons that prevented detecting of the 
radio signals. 

Considering Eq. 1 1 in Iacolina et al. (2009), we estimated 
the probability that the beams of the sources do not intersect our 
line of sight. Assuming a typical duty cycle of 15% for each 
source, the probability that the beams of all the sources (includ- 
ing XTE J0929-314) are missing our line of sight is ~19%. In 
order to exclude this geometric bias we would have to analyse 
the whole known sample of twelve AMXPs, and the probability 
of missing all the beams of the whole sample will drop to about 
0.1%. 

The second possible reason that might have prevented the 
detection of a radio signal could be the low luminosity of the 
sources. As the true luminosity of pulsars cannot be measured 
because of the unknown beaming fraction, a 'pseudo luminos- 
ity', L, is defined as the observed flux density, S, multiplied by 
the square of the pulsar distance, d (Taylor & Manchester 1977): 
L = S x d 2 . 

In Fig. 2 the logarithmic distribution of the pseudoluminos- 
ity at 1.4 GHz of the sample of 46 known galactic field MSPs 
is shown 1 . The vertical red lines indicate the lower values of 
the upper limits on L of XTE J1751-305, XTE J1814-338, 



1 Data taken from the ATNF pulsar catalogue - http : //www . atnf . 
csiro.au/research/pulsar/psrcat/; Manchester et al. (2005). 
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XTE J1814-33B 
XTE J1751-305 



Log [Sxd s (mjy kpc 2 )] 



Fig. 2: Pseudoluminosity distribution of a sample of 46 MSPs in the 
galactic field. The red lines indicate the upper limits related to the four 
analysed sources. 



SAX J1808.4-3658, and XTE J0929-314 (for the observations 
at 8.5 GHz and 5 GHz), scaled at 1.4 GHz, considering a depen- 
dence on the frequency S(v) oc y~'\ with an index a — 1.7, and 
for the distances indicated in Table 3. 

These limits determine the probability that the true pseudo 
luminosity of each source is too faint for detection in our search. 
We calculated that this probability is about 90% for XTE J1814- 
338 and XTE J1751-305 and -80% for SAX J1808.4-3658 and 
for XTE J0929-314. The combined probability of the non de- 
tection due to the luminosity threshold of our survey is ~50%. 
This percentage is not negligible and can be reduced by a deeper 
search and/or by a larger sample. 

Expanding the sample to all twelve known AMXPs, with a 
probability equal to about 80% for each one, the combined prob- 
ability would be ~10%, not enough for a safe detection. For a 
combined probability less than ~0.1%, we also have to perform 
a deeper search. The upper limit in pseudo-luminosity for each 
source having such a probability is L m 2 mjy kpc 2 . The limit 
in flux for an average distance of 7 kpc becomes S ~ 0.04 mjy 
at 1.4 GHz which, scaled at 4.7 GHz (as, for example, for the 
observation of SAX J1808.4-3658), becomes S ~ 0.003 mjy. 
Such a limit can be reached by performing observations using 
telescopes with a larger bandwidth and higher instantaneous sen- 
sitivity. In fact, a 4.7 GHz observation of SAX J1808.4-3658, 
performed using a 2 GHz bandwidth and 2.01 K Jy~' gain ob- 
tainable at the Green Bank radio telescope (GBT), would have 
reached a flux density limit of a;0.003 mjy. 



Possible mechanisms that might have hampered the observa- 
tion of the pulsed signal could concern the luminosity of the four 
analysed sources, lower than our limit of detection, resulting in a 
~50% combined probability of non detection, or the anisotropic 
nature of the pulsar emission, with a probability of ~19% that 
the beam of all the four sources does not intersect our line of 
sight. 
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5. Conclusions 

We performed a search for millisecond pulsations in three 
AMXPs, XTE J1814-338, XTE J1751-305, and SAX J1808.4- 
3658, in their quiescence phases at high radio frequencies. 
Discussion was done by including the previous results obtained 
by Iacolina et al. (2009) for XTE J0929-3 14. Except for the case 
of XTE J175 1-305, for which further investigations are needed, 
no pulsations with the expected periodicity have been detected 
in the analysed data. The flux density upper limits determined by 
our search (including XTE J0929-314) are indicated in Table 2. 
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Table 3: Parameters for the source, the observation and data analyses for XTE J1814-338, XTE J1751-305 and SAX J1808.4-3658. 



Pulsar names 


XTE J1814-338 


XTE J 175 1-305 


SAX J 1808. 4-365 8 


Source parameters 


Right ascension (J2000.0) 


18 h 13 m 39?04 


17 h 51 m 13f49 


18 h 08 m 27f54 


Declination (J2000.0) 


-33°46'22'.'3 


-30°37'23'.'4 


-36°58'44'.'3 


Orbital period, P mb (s) 


15388.7229(2) 


2545.342(2) 


7249.156980(4) 


Projected semi-major axis, a sin i (lt-ms) 


390.633(9) 


10.125(5) 


62.812(2) 


Eccentricity, e 


< 2.4 x 10~ 5 


< 1.3 x 10~ 3 


< 1.2 x 10~ 4 


Spin period, P s (s) 


0.0031811056698(1) 


0.002297172972(2) 


0.00249391975978(6) 


Mean spin period derivative, P$ (s s~') 


6.7(7) x 10~ 19 


-1.95(10) x 10~ 18 


4.7(9) x 10~ 19 


Ascending node passage, T (MJD) 


52797.8101689(9) 


52368.0129023(4) 


52499.9602472(9) 


Distance, d (kpc) 


8 


6.7 


3.5 


Nominal DM (pc cm 3 ) 


-200 


-400 


-100 


References 


1 


2 


3,4 


Observation parameters 




j st 2 nd 




Central radio frequency, y b s (GHz) 


6.4105 8.4535 


8.4535 8.4535 


4.7495 6.3515 


Data series time, At (s) 


25053 10669 


17145 22864 


35685 35051 


Sampling time, St ssm (/us) 


92 50 


60 100 


125 125 


Samples number, N s 


2 28 2 27 


2 28 2 27 


2 28 2 28 


Nominal gain , G (K Jy" 1 ) 


0.46 0.59 


0.59 


0.46 


System temperature" , r sys (K) 


50 25 


25 


50 


Number of 3 MHz channels, N c 


192 


192 


192 


Data analysis parameters 


Number of DMs 


186 151 


195 127 


135 57 


DM range (pc cm 3 ) 


-50-1000 


- 100-1700 


- 20-400 


Local DM (pc cm 3 ) 


-20 


-260 


-230 


Number of /> or b steps 


1 


19 


1 


Number of folding steps 


1 


11 


1 


Other parameters 


Sky temperature, T sky (K) 


0.06 0.03 


0.09 


0.13 0.06 


Dispersion broadening, <5? DM (10~ 8 s) 


9.5 x DM 4.2 x DM 


4.2 x DM 


23 x DM 9.7 x DM 


Scattering broadening, <5/ scatt (s) 


9.3 x IQ- 9 3.2 x 10-" 


4.4 x 10~ 7 


9.6x10-* 3xl0- 9 



Notes. (a) Values derived from the Parkes website: http : //www.parkes . atnf . csiro . au/observing/documentation/. 
References. (1) Papitto et al. (2007); (2) Papitto et al. (2008); (3) Hartman et al. (2009); (4) Burderi et al. (2006). 



